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Chirally odd generalized parton distributions (GPDs) and the Boer-Mulders function provide 
valuable information about spin-orbit correlations for quarks in nucleons and other hadrons. We 
compare results for the relevant GPD from a variety of phenomenological models as well as 
recent lattice results. We find that £^ > for nucleons as well as the pion and for both u and 
d quark. As a result, the corresponding Boer Mulders functions are all expected to be negative. 
The sign of arises from the relative sign between the upper and lower Dirac components for the 
quark wave functions. 



I. INTRODUCTION 



Recent experiments by the Hermes collaboration have demonstrated the existence of significant single-spin asym- 
metries (SSAs) in semi-inclusive deep- inelastic scattering (SIDIS) The Sivers function [5] measured in these 
experiments indicates a negative correlation between the transverse proton polarization and the transverse momen- 
tum for ejected u quarks, while the same correlation for d quarks turned out to be positive (Note that we adhere to 
the Trento convention fsj). 

Meanwhile, significant advances have been made in the theoretical understanding of SSAs For example, simple 
model calculations have illustrated very clearly that SSAs do not need to vanish in the Bjorken limit [5*1 . However, it 
has also become clear that several ingredients need to conspire in order to produce SSAs ^6(]: essentially what is needed 
is an interference between final state interaction phases for quarks originating from different partial waves in the target 
wave function. Clearly, such a complex interplay between different aspects of hadron structure not only makes SSAs 
interesting, but it also presents a challange for theorists to produce predictions that do not depend sensitively on 
model parameters. Indeed, due to the complexity of these caclulations, it even took a while until agreement was 
reached in observables as simple as the sign of the Sivers function in the scalar diquark model. 

Significant clarification regarding the sign issue was accomplished once it was noticed that the sign of these SSAs 
can be intuitively related to the sign of the deformations of certain spin-dependent impact parameter dependent 
parton distributions 0, Q: the final state interactions (FSIs) are expected to be attractive, on average, and therefore 
impact parameter dependant parton distributions (IPDs), which are purely intrinsic properties of hadrons, can be 
related to SSAs, which also invole FSIs. In the case of the Sivers effect, the relevant IPD is the Fourier transform 
of the Generalized Parton Distribution (GPD) E"^ {x , , — A.\) that describes the ^-decomposition of the Pauli form 
factor F2(— A^^). Since the latter is approximately known for q = u,d, it was possible to predict the sign of the u 
and d quark Sivers functions [§, [l^, [ij, [l2| which were subsequently confirmed by the Hermes data ^ . 

Even when the target nucleon is unpolarized, the momentum distribution of its quarks ejected in SIDIS can already 
exhibit a left-right asymmetry of the transverse quark momentum relative to their own transverse spin Sq [13] 



fqyp{x,kT) = i 



2^ ,.2^(PxkT)-S, 



(1) 



P is a unit vector opposite to the direction of the virtual photon momentum and M is the target mass. This 
asymmetry can be measured experimentally by tagging the transverse polarization of quarks produced in SIDIS using 
the Collins effect Hf]. Alternatively, one can also study the cos 20 asymmetry in the unpolarized Drell-Yan process, 
where </> is the azimutal angle of the ^~ plane about the virtual photon axis w.r.t. the incident proton. Like the 
Sivers function, the "Boer-Mulders function" h^'^{x, k'^) requires a combination of orbital angular momentum with a 
difference between the final state interaction from different parts of the initial quark orbit. Since the Boer-Mulders 
function specifies the correlation between the transverse momentum asymmetry of the struck quark and its spin, it 
provides important information about the correlation between quark orbital angular momentum and spin. Similar 
to the case of the Sivers function, one expects that the sign of h^'^ can be related to the distribution of transversely 
polarized quarks in impact parameter space [15l |. which can be obtained from the Fourier transform of chirally odd 
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GPDs 
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+ + 2Ht)u + Erix, C t)u u 



The linear combination Et = Et + 2Ht multiplying the last term in ([3]) describes the transverse deformation of the 
distribution of transversely polarized quarks in an unpolarized nucleon JJJJ . In ref [l5l | it has been suggested that h^'^ 
and E!^ have opposite signs. Therefore, knowledge of E!^ make qualitative predictions for the Boer-Mulders function 
possible. However, not even the sign of E!^ is known experimentally. The main purpose of this paper is to develop 
some intuition about the expected sign for Et by analyzing this observable in a variety of hadron models. 

II. IMPACT PARAMETER DEPENDENT PDFS IN THE BAG MODEL 

The well known MIT bag model offers a familiar venue to begin our discussion of the sign of SSAs in the nucleon. 
Bag model wavefunctions have the form 



N 



where Xm is a Pauli spinor, and the jn are the spherical Bessel functions and iV is a normalization factor. In order to 
slightly broaden the scope of our discussion, note that ji related to jo by a derivative. In fact, in the bag model, this 
relation between the upper and lower component is a consequence of the (free) Dirac equation. With this in mind, 
let us instead focus on a more general Dirac spinor, 

^rn = ( .^^"a. ) ' (5) 



-gi^ ■ x)Xr, 

where / is a monotonically decreasing radial function, and g is the derivative of /. 

The impact parameter dependent parton distributions that we would like to evaluate are of the form 

Fr{x,h^)=Af-^ j ^e"f+^" (p+,Oi|g(0,bi)rg(z-,bi)|p+,Oi). (6) 

Here |p+, 0_l) is an eigenstate of light-cone momentum p+ as well as of the transverse center of momentum operator 
Rj^ (for a definition of the latter see Refs. 0,[l3). 

From the bag model wave functions one can easily evaluate position space densities of quark bilinears, but evaluating 
parton distributions requires introducing additional recipies in order to deal with light-like correlation functions. This 
issue can be entirely avoided by focusing on the lowest moment of GPDs, yielding a (local) density that only depends 
on r^. 

JdxFr{x,h^) = (2p+AA)-i(p+,Ox|g(0,bi)r<z(0,bi)|p+,Oi). (7) 

Instead of evaluating the matrix element of the operator in Eq. ([7]) at the origin between plane wave states, one may 
equivalently localize also the longitudinal position of the state and instead integrate the operator over longitudinal 
position, yielding 

J dxFr{x,h±) = const. J dx^ {6\ q{x^ ,h±)rq{x^ ,h±) \b) , (8) 

where |o) denotes a nucleon state localozed at the origin (in all three space directions). Note that, as the bag model 
is not boost invariant, there is a certain arbitrariness in the extraction of Et from bag model wave functions. In 
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particular, the reference point for the impact parameter bj^, which should be the transverse center of longitudinal 
momentum, is here taken to be simply the center of the spherical bag. However, since we are only concerned about 
the overall sign of Et, such issues should not be significant. 

Quarks with transverse polarization s are projected out by the operator [7^^ — sHa~^^^5j q [l6j and therefore the 
vector field representing the transverse quark polarization density is given by —iqa^^^f^q. We thus consider impact 
parameter dependent PDFs with F = —ia^^j^ which are related to the (Fourier transforms of the) chirally odd GPDs 
Et, Ht and Ht [H 

= -^''f^^^T + (ht 4^A,7iT) + {2b'V - b'S^^) S^^K (9) 

where we used script letters to denote the Fourier transforms of these GPDs and is the spin of the target. Only 
the term involving £t contributes for an unpolarized target (averaging over transverse target polarization), which is 
why it is only the GPD E^ that is expected to be (qualitatively) related to the BM function. In the bag model, we 
extract this term by considering the density corresponding to F = —ia^^^5 and summing over the target spin. For a 
single quark state this procedure yields 

5^(P5™|^'(a;^bi)^a+^5*(x^bi)|P5™) = -^Y.^f + d'V-^ + ^^9^''^^^ ~ ^9'b\{b^ ■ s,^,) (10) 

where Sm is the spin vector corresponding to the pauli spinor Xm- The first and last terms of (jlOp do not survive the 
sum over 'target' polarizations. The asymmetry is given entirely by the middle term, which is an interference between 
the uppper and lower components of ([5]). For the lowest moment of £t, we find 
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KT ^ I dxET{x,{),Q) ^ j dxd?h^£T ^ I drr^fg. (11) 
J J 3\/2'K Jo 

The right hand side of pip is always positive because / and g are non-negative functions for r less than the bag 
radius, implying that > 0. 

In the bag model, the correlation between quark spin and quark orbital motion is the same, regardless of the 
orientation of jz- AH quark spin orientations thus contribute coherently to and in the case of d quarks, £^ is equal 
to for a single quark, while for u quarks it is twice as large. In fact, for any model where the quarks are confined 
by some mean field potential one finds that all quark orbits give the same contribution to and thus £!^ is equal 
to £^ for a single quark orbit, multiplied by the number of quarks of flavor q. In particular, in the large Nc limit, 
where = + 1 ^ 00, the lowest x moment of is the same for u and d quark and both are of order 0{Nc)- 
Since the support of GPDs shrinks to x ^ 0{1/Nc), this implies that E^{x,^,t) = E^{x,^,t) = 0{N^) 

In order to visualise the transverse spin - impact parameter correlation in the bag model, the vector field 

dx^fge'^b^ (12) 

representing the lowest moment of the transversity density in an unpolarized target has been plotted in Fig. [l]for bag 
model wave functions / = joif), and g = ji{f)- In the bag model, we thus optain a counter-clockwise polarization for 
impact parameter dependent quark distributions. Upon invoking the "chromodynamic lensing" mechanism (attractive 
FSI), and using the Trento convention 0], this implies a negative Boer-Mulders function h^. This result is also 
consistent with a direct calculation of for the bag model in Ref. [l^. The most important question here is: how 
general is this result regarding the sign of h^, i.e. how much does it hinge on specific features of the bag model. 

In order to address this issue, let us consider the general case of a particle confined by a combination of scalar 
potential (i.e. a mass term m{r) that depends on the radius) as well as a (zero component of a) vector potential V{r). 
The bag model corresponds to the limiting case where the vector potential vanishes and the scalar potential has the 
shape of an infinite square well. For a free quark, the lower component (pi is related to the upper component 4>u, via 
the free Dirac equation 

-a-p(j)u. (13) 



" E- 

In the presence of an external potential this relation changes only slightly 

1 ^ ^, 



E + m{r) - V{r) 



(14) 
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FIG. 1; Lowest moment of the impact parameter dependent transversity distribution for an unpolarized target in the MIT bag 
model (I12|l . The 'outside' of the spherical bag corresponds to the regions without arrows. 



The vector part of the confinement potential cannot exceed the scalar part — otherwise one encounters the Klein 
paradox. Therefore, E + rn — V should in general be positive As a result the sign in this relation is unchanged by 
the presence of the potentials m{r) and V{r). In particular, since one would expect that the upper component /(r) 
for the ground state is a monotonically falling positive function (its overall phase has been taken positive here), one 
expects g{r) to be negative, regardless of the details of the potentials. Note that while E + m — V can be negative 
at isolated points without giving rise to the Klein paradox, e.g. very close to the origin for an electron moving in a 
Coulomb potential, this does not invalidate our main point regarding the average value of the product /(r) • g(r) or 
/(r) • g{r)r relevant for the lowest moments of Et at zero momentum transfer. 

This discussion illustrates that the sign of the spin-orbit correlation (counterclockwise) described by Eq. (jlOp should 
be the same for the ground state of all confining potential models. 



III. DIQUARK MODELS 



One of the main disadvantages of the bag model is its lack of boost invariance, which results in certain ambiguities 
in the parton interpretation of bag model wave functions. In diquark models the distribution of quarks in a 'nucleon' 
is obtained by allowing it to split into a quark and a spectator 'diquark'. As the wave function is usually calculated 
perturbatively using a point-like vertex, boost invariance and electromagnetic gauge invariance are straightforward to 
maintain. The same applies to (perturbatively calculated) QCD final state interactions between the quark and the 
diquark. Because of the latter features, it is thus not surprising that the scalar diquark model has provided the first 
clear and convincing example for the fact that SSAs can survive in the Bjorken limit [5,]. 

The phenomenological motivation for diquark models is the idea that, as long as the momentum transfer on the 
spectators is not very large, one may replace them by a single, point-like degree of freedom — the diquark. The 
quark-quark interaction in the nucleon is more attractive in the scalar channel with quark spins anti-alligned. This 
further motivates the assumption that the spectator diquark carries scalar quantum numbers. In this work, we will 
mostly focus on the scalar diquark model only, although more sophisticated diquark models exist. 

The scalar diquark model wave function, has only two spin components corresponding to positive (t) and negative 
quark helicity (J,). Angular momentum conservation implies that the wave function components for 'nucleon' states 
with positive (ff) helicity take on the form 



= /(a:,ki) 

= (r +zF)g(x,ki) 



(15) 
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Time reversal implies for the negative helicity states (-IJ-) 

i^f{x,-k^) = {-k- + tky)g{x,kl) (16) 
^lix,k^) = f{x,kl). 

For transversely polarized quarks (lower index) in a transversely polarized target (upper index) this implies 

^XUx,k^) - /(x,ki)+zFg(x,ki) (17) 
i^ll{x,k^) = -rg(x,kl), (18) 



i.e. in impact parameter space 



^+|(x,Ci) ^ J d2k^e*---^+|(x,k^) = f{x,cl) + ^~9ix,cl) (19) 



J d2k^e*---V^I(^,k^)=^^ff(a;,ci) 



As a result of the interference between the helicity non-flip / and hclicity-flip component dyg, the distribution of quarks 



with the same transverse polarization as the 'nucleon' 



exhibits a transverse asymmetry proportional to 



fdycj. The wave function component where the transverse quark spin is opposite to the nucleon spin does not exhibit 
such an asymmetry. Intuitively, one can understand the origin of the asymmetry in the scalar diquark model as follows. 
When the 'nucleon' splits into a quark and a diquark, the orbital angular momentum of the two particle system can 
be either ^ = or ^ = 1 (Note that this is not in conflict with parity invariance as upper and lower components of 
quark spinors transform differently under parity and leading twist quark densities result from superpositions of upper 
and lower components). Higher orbital angular momenta are not possible as the total angular momentum (0AM plus 
spin of the quark) must still be ^. For quarks that have the same (transverse) spin as the 'nucleon', the orbital wave 
function can have either / = or Z = 1 and Ix = 0. The asymmetry results from the interference between these two 
wave function components very much in the same way as hybridization in chemistry results in asymmetric electron 
orbitals. In contradistinction, the wave function component where the quark spin is opposite to the 'nucleon' spin 
is pure 1 = 1, i.e. there cannot be any interference between I = and I = 1 and hence no asymmetry. This very 
simple observation also 'explains' why the BM and Sivers functions are identical in the scalar diquark model (for an 
explicit calculation of both functions in this model see Ref. [l^): a left-right asymmetry only occurs in the wave 
function component where the transverse spin of the quark and the transverse spin of the 'nucleon' are the same. 
Thus the correlation between the transverse position space asymmetry with the quark spin is the same that it is 
with the nucleon spin. Since (at least in diquark models) the transverse momentum space asymmetry in SIDIS arises 
entirely from the initial transverse position space asymmetry, the BM function (representing the correlation between 
transverse momentum and spin of the struck quark) and the Sivers function (representing the correlation between 
transverse quark momentum and 'nucleon' spin) must be the same. While it has been known that those two T-odd 
functions are the same in the scalar diquark model, the intuitive explanation provided here is new. 

For models containing axial vector diquarks the situation is mode complicated, as there can also be interference 
between wave function components where the transverse quark spin is opposite to the 'nucleon' spin (and the vector 
diquark spin is along the 'nucleon' spin direction). As a result, the correlation between quark spin and its transverse 
position is no longer identical to the correlation between 'nucleon' spin and quark transverse position. Hence the 
Sivers and BM functions in vector diquark models do not have to be identical. 

In the scalar diquark model 



u+y( 



m Id 
M+-+ -— 
X X del' 



= -^^^(-l) (20) 



where is the distance between the active quark and the spectator and 0(cj^) is monotonically decreasing with c'^. 
This result is of the above form p^ . with / = [M + ^] and g = ^0, i.e. quarks polarized in the +x direction 
are shifted towards negative y. Again, we find a counter-clockwise pattern for the transverse quark polarization 
corresponding to a positive GPD E!^. A qualitatively similar result is obtained for axial vactor diquark models. 

It is easy to understand why we obtain the same sign as in the bag model. The quark interacts with the diquark (to 
form the nucleon) only when the two are at the same point, i.e. in impact parameter space the interaction is a contact 
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interaction. For nonzero impact parameter, the quark is thus a free particle and as such subject to the free Dirac 
equation relating upper and lower component. Just as was the case in the bag model, this determines the relevant 
phase which in turn determines the sign of Ej^. This observation has far reaching consequences for other (e.g. axial 
vector) diquark models. What all diquark models have in common is the fact that quark and diquark interact only 
when they merge to become the nucleon, i.e. the quark-diquark interaction is local. Both for scalar and axial vector 
diquarks, the the quark-diquark state is mainly an s state and therefore the same general statements regarding the 
sign of i?^ apply to both the axial vector and scalar case. 

IV. RELATIVISTIC CONSTITUENT QUARK MODELS 

In relativistic constituent quark models for hadron structure, one starts from nonrelativistic forms for the quark 
wave function. Nontrivial spin structure is then generated by boosting these wave functions to the infinite momentum 
frame. In this procedure, using non-interacting boost operator, the net result is that instant form spin eigenstates 
\k^,'k±;±)^ are replaced by light-front helicity eigenstates |fc+,kj^;±)^ using what is often referred to as a Melosh 
transformation 

|fc3,kx;+), = u;[(fc++m)|fc+,ki;+)^-fc^|fc+,ki;-)^] 

|fc^kx;"), - u;[(fc++m)|fc+,kj.;-)^-f fc^|fc+,k_L;+)^] (21) 

where w = l/-\/2fc+(m + fc*^), k^/^ = k^zLik'^, and A:+ = k'^ + k^. For quarks with nonzero transverse momentum, the 
boost to the infinite momentum frame (or equivalently the Melosh transformation) naturally introduces a corelation 
between quark spin and quark orbital angular momentum and a non-trivial spin structure of the nucleon state — 
even if the original (instant form) wave function only contained s-wave components and SU(6) wave functions. 

Such models have been used to estimate E!^ in Ref. ^21j|. Studying a wide range of model parameters, it was found 
that in the range between 1.98 and 3.60, and Kj_ in the range between 1.17 and 2.36. The sign of these results is 
in agreement with the results from the other models. Furthermore, similar to the bag model, one finds Ku ~ 2^^. It 
is possible to understand the origin of the sign for Et in this model. The relativistic constituent models start from 
pure s wave wavef unctions. Any p wave component is generated by the Melsoh transformation — which employs the 
free (noninteracting) boost operator. It is thus not surprising that the phase relation between the s and p wave is 
thus the same as for all the other models studied here. 

V. Et in the PION 

Since the pion has no spin, its impact parameter dependent PDFs are much simpler than in the case of the nucleon. 
The distribution of quarks with spin s' in impact parameter space reads 

i [F + s^F^] ^H{x,h') + s\'^b^^^ET{^,h% (22) 

where 7i(a;,b^) and £j'(a;,b^) are again the fourier transforms of the GPDs H{x,0,t) and ET{x,0,t) respectively 
whose definition is particularly simple 

/ (V|g-(-lz)7+9(iz)gk) U.=o,.=o = Hix,^,t) (23) 

/ (^'1 q{-lz)a+^j,q{^z) \n) U+^o,.=o = ^^t(x, '''' p^"^' • (24) 

Here A is some hadronic mass scale, which needs to be included in the definition if one wants ET{x,£,,t) to be 
dimensionless (cf. Eq. [3]). For the nucleon, the natural choice is the nucleon mass, but chosing A = m,r for the pion 
would perhaps not be a wise choice, as this would unnecessarily complicate the discussion of the chiral limit. Since 
we are again mainly interested in the sign of Et, we leave this issue open, but perhaps A ~ 4:7: f^^ or A = A/jv would 
be more useful choices for the pion case. 

Except for a slight change in the bag radius, the quark wave functions in the bag model are the same for pions and 
nucleons. Therefore, apart from a slight rescaling due to the different bag radii, EJj^ in a 7r+ is the same as -^E^ or 
E^ in the proton. Here the factor i accounts for the fact that there are twice as many u quarks in a proton than in 
a 77+. 
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As an alternative model for GPDs in the pion we also considered the Nambu-Jona-Lasino (NJL) model. In this 
model, the quark and the anti-quark interact via a 75 contact interaction in the s-channcl, and form factors or GPDs 
are obtained by evaluating momentum intergrals containing quark traces of the form 



tr [75 (fc + f + M)r(fc - f + M)75(fc -p + M)] 



{k. 



M2 



[{k-py-M^] (fc 



AP 



(25) 



Here M is the (constituent) quark mass and T depends on the current under consideration. For example, T — 7+ is 
used to calculate PDFs and the GPD H, yielding for the trace 



tr[r = 7+] = 4 {p+(A/2 _ fc2) + fc+[Af2 - (p - fc)2 + p2j 



(26) 



Evaluating the expression when the spectator is on-shell [{p — fc)2 = M2], which arises from complex countour 
intagration for k~ , and dropping all terms proportional to rn^ this yields 



tr[T = 7+] ^ 4p" 



M2 



.Af2 



fc+ 



1-x 



(27) 



Together with a factor in the energy denominators, this gives rise to the well known result of a constant PDF in 
the NJL model (for to^ — 0). 

In the calculation of AiEx, it is F = 7'^7' which enters the above trace. For simplicity, we focus here on Et 
for A —> which determines the average spin dipole moment kj_. Since the denominator is even in A, only the 
A-dependence in the numerator matters in this limit, and we thus consider only the relevant numerator 



tr[T = 7+7'] = tr 



(fc + ^ + M)7+7'(fc -^+M){p^k + M) 



^A,Mp+, 



= [[(fc + Af)7+ + 7+(fc + M)]ip - fc + Af )] 

(28) 



yielding again > 0. A more detailed analysis yields Ex oc (1 — a;), for TOtt 0, where the constant of proportionality 
depends on the cutoff in the loop integrals. 

Enlightened by our discussion about the sign of Et in the nuclcon, it is easy to uderstand the result in the NJL 
model: the interaction between quark and antiquark in this model is also a contact interaction, i.e. the quark and 
antiquark follow the free Dirac equation, except when they sit on top of each other. Therefore, by the same reasoning 
as in the bag model, Et should be positive, as the explicit calculation confirmed. Lattice calculations yield the same 
sign (2^ . 



VI. SUMMARY AND CONCLUSIONS 



We have studied the chirally odd GPD Et for both u and d quarks in a variety of models (Bag model, potential 
models, diquark model, NJL-model, relativistic constituent quark models). All models provide the same sign Et > 0. 
This sign also agrees with recent results from lattice QCD calculations. While the physical origin of the sign is 
obscure in lattice QCD calculations, it is much more transparent in the model calculations. As the quark helicity flips 
in the matrix elements contributing to Et, the nucleon helicity does not. Hence the quark orbital angular momentum 
between the initial and final state must differ by one unit and therefore Et arises from the interference between wave 
function components that differ by one unit of orbital angular momentum. In many models for the structure of ground 
state hadrons only s and p wave components are present and higher 0AM is either completely absent or negligible. 
Moreover, in these models, the p-wave component arises primarily from the 'lower' component of the hadron wave 
function, which is related to the upper component via the Dirac equation, which thus determines the relative phase 
between the upper and lower components. While the presence of interactions tends to modify the relation between 
the upper and lower components quantitatively, it does not change the phase relation and all models yield the same 
phase relation between upper and lower component that also holds for free quarks. Lattice results confirm this result. 
Of course, one needs to be cautious as present day lattice calculations are still performed at moderately high quark 
masses. The pion masses used in Ref. 23 1 range from « 1 GeV^ to « 0.2 GeV^. and the fact that the phase 
relation between upper and lower components is the same as for the free Dirac equation could in principle be an 
artifact of the large quark mass. However, the results from Ref. [1^ are also rather stable against variations of the 
quark mass, and the sign of E!^ does not change upon extrapolation to the physical quark masses, yielding w 2.93 
and « 1.90 
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Such a unanimous agreement between such a variety of theoretical approaches to an observable that has never 
been measured is probably almost unprecedented. While some of the model predictions for the sign of existed 
before, what is new in this work is the intuitive understanding about the sign of Ej. as resulting from the relation 
between upper and lower components in the (free) Dirac equation. Despite all the complexities of QCD, lattice 
calculations yield results that are qualitatively consistent with quarks in an s state that have a p-w&ve in the lower 
(Dirac) component of the wave function and with a relative phase obtained from the free Dirac equation. While this 
agreement may be a pure accident, it is still very tempting to interpret this result as an indication that a nontrivial 
fraction of the orbital angular momentum in the nucleon wave function is merely resulting from the p-wave admixture 
in the lower component of quark spinors in a bound state. 

Acknowledgements: This work has been partially supported by the DOE under grant number DE-FG03- 
95ER40965. 
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